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Abstract 

The self-consistent proton-neutron quasiparticle random phase approximation approach is em¬ 
ployed to calculate /3-decay half-lives of neutron-rich even-even nuclei with 8 Z ^ 30. A newly 
proposed nonlinear point-coupling effective interaction PC-PK1 is used in the calculations. It 
is found that the isoscalar proton-neutron pairing interaction can significantly reduce /3-decay 
half-lives. With an isospin-dependent isoscalar proton-neutron pairing strength, our results well 
reproduce the experimental /3-decay half-lives, although the pairing strength is not adjusted using 
the half-lives calculated in this study. 
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The nuclear /3 decay is an important property of nuclei, which can play an important role 
in nuclear physics and astrophysics. In nuclear physics, the investigation on /3 decay can 
provide information on spin and isospin properties of nuclear effective interaction as well as 


nuclear properties such as masses 


Cl, shapes [2:], 


and energy levels |3| • In astrophysics, about 


half of the nuclei heavier than Fe are synthesized by the rapid neutron-capture process (r 
process). The /3-decay rates of nuclei on the r-process path set the time scale of the r-process, 
hence determine the production of heavy elements in the universe C l|- 

With the development of radioactive ion beam facilities, the measurement of nuclear 
/3-decay half-lives has achieved great progress in recent years d • For example, the measure¬ 
ments of the /3-decay half-lives of neutron-rich nuclei in the vicinity of N = 28 and N = 50 


shell closures 


11}. On the theoretical side, apart from the macroscopic gross theory 12], 


two different microscopic approaches have been employed in the theoretical predictions of 
nuclear /3-decay rates, i.e., the shell-model and the proton-neutron quasiparticle random 
phase approximation (QRPA). 

The nuclear shell model can take into account the detailed structure of the /3-strength 
function [6], while the QRPA approach provides a systematic description of /3-decay prop¬ 
erties of heavy nuclei. The early phenomenological QRPA approaches for the /3-decay half- 
lives rely on separably effective nucleon-nucleon interactions, such as the QRPA approach 


based on the deformed Nilsson+BCS formalism 


13] and on the finite-range droplet model 
Recently, the self-consistent 


(FRDM) with a folded Yukawa single-particle potential [Q. 

QRPA approach has received more attention, which is generally believed to possess better 
extrapolation ability than the phenomenological QRPA. 

The self-consistent QRPA approach has been employed in the calculation of nuclear 
/3-decay half-lives based on the extended Thomas-Fermi plus Strutinsky integral (ETFSI) 


model 


15}, the Skyrme-Hartree-Fock-Bogoliubov (SHFB) model 16], or the density func¬ 


tional of Fayans 17]. On the other hand, the relativistic model has also received wide 
attention due to many successes in describing lots of nuclear phenomena [18|, Il9| and suc¬ 
cessful applications in astrophysics 20-122], In the relativistic framework, there are two 
widely used models, i.e., the finite-range meson-exchange and zero-range point-coupling 
models. In the meson-exchange model, the QRPA has been developed based on the rela¬ 


tivistic Hartree-Bogoliubov (RHB) approach [23] and employed in the calcu 


half-lives of neutron-rich nuclei in the N ~ 50 and N m 82 regions 24J, [25]. Moreover, the 


ation of /3-decay 
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self-consistent relativistic QRPA was also developed based on the relativistic Hartree-Fock- 


Bogoliubov (RHFB) theory (26|, 27] and it has been used to calculate /3-decay half-lives of 


neutron-rich even-even nuclei with 20 ^ Z ^50. With an isospin-dependent T = 0 proton- 
neutron pairing interaction, the known half-lives are well reproduced in the whole region of 
20 ^ Z ^ 50 

In the point-coupling approach, the proton-neutron QRPA has been derived based on the 
effective Lagrangian with density-dependent couplings determined by chiral pion-nucleon dy¬ 
namics 28]. In addition, the nonlinear point-coupling effective interactions have also received 


wide attention due to its successful description of lots of nuclear phenomena 


29]. Recently, 


a new nonlinear point-coupling effective interaction, i.e., PC-PK1, was proposed, which well 
reproduces the properties of infinite nuclear matter and finite nuclei including the ground- 
state and low-lying excited states 


30 


33}. Based on the nonlinear point-coupling effective 


interaction PC-PK1, the self-consistent QRPA has been formulated and was employed to 
calculate the /3 + /EC-decay half-lives of neutron-deficient Ar, Ca, Ti, Fe, Ni, Zn, Cd, and 
Sn isotopes |34|. The isospin-dependent pairing strength proposed in Ref. 27] predicts a 
constant values of Vo for these neutron-deficient nuclei and the calculations found that their 
half-lives are indeed reproduced well by an universal pairing strength 34], 


In this work, we will calculate the /3-decay half-lives of neutron-rich nuclei with the newly 
proposed effective interaction PC-PK1. The isospin-dependent pairing strength is adopted 
but the parameters are directly taken from those in Refs, 
with the known half-lives of calculated nuclei in this work. 


27], l34|, i.e., not again adjusted 


The relativistic proton-neutron QRPA has been formulated in Ref. [23| using the canonical 
single-nucleon basis of the RHB model for the meson-exchange effective interaction. Based 
on point-coupling relativistic Hartree-Bogoliubov theory, the self-consistent proton-neutron 
QRPA has also been established recently 28], |34], Similar to Ref. 34], isovector-vector 


interaction, direct one pion interaction and the corresponding zero-range counter term are 
included in the QRPA particle-hole (p-h) residual interaction. The specific expressions for 
these terms can be found in Ref. 34]. 

For the particle-particle (p-p) residual interaction, the isovector (T = 1) and isoscalar 
(T = 0) proton-neutron pairing interactions are included in the QRPA calculations. In the 
T — 1 p-p channel, we employ a phenomenological pairing interaction, i.e., the pairing part 
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of the Gogny force 


V T = i( 1, 2) = e- [(ri ” r2)/w]2 

1 = 1,2 

(Wi + B % P° - H t P T - M i P ff P T ), 


( 1 ) 


with the parameter set D1S for /q, H+ B^ H ll and Mi 35]. For the T — 0 proton-neutron 


pairing interaction in 


Refs. 


16 


24 


25 


m 


34J: 


he QRPA calculation, we employ a similar interaction as in the 


R t=0 (1,2) = —R 0 ^,, 


,-[(ri-r 2 )//ij] : 


n 


3 =1 


S=1,T=0 


( 2 ) 


with pi = 1.2 fin, /i 2 = 0.7 fin, pi = 1, g 2 = —2. The operator J][ s=1 r=o P r °j ec ts onto states 
with S — 1 and T = 0. 


The strength parameter V 0 in Eq. ([2]) is usual 
of selected nucleus in each isotopic chain 


24 


y determined by fitting to known half-lives 


25| . However, very different values are found 


for different isotopic chains, which limits the prediction power of the model. For improving 
this dilemma, an isospin-dependent pairing strength Vo was proposed in Ref. [27|, i.e., 

V D 


Vn = VV 


(3) 


1 _|_ e a+b{N-Z) ’ 

where V L = 134.0 MeV, V D = 121.1 MeV, a = 8.5, and b = —0.4 are obtained by fitting to 
the known half-lives of neutron-rich nuclei with the calculations of RHFB+QRPA approach. 
This isospin-dependent pairing strength predicts a constant values of Vg for nuclei with N — 
Z < 5, which is just the case for the neutron-deficient nuclei in the region 20 ^ Z ^ 50. For 
further testing the validity of this expression, the half-lives of neutron-deficient Ar, Ca, Ti, 
Fe, Ni, Zn, Cd, and Sn isotopes have been calculated with RHB+QRPA approach and PC- 


PK1 effective interaction 


341 ]. It is found that their half-lives can be well reproduced by an 


universal pairing strength Vo = 175.0 MeV, which indicates the validity of this expression of 
Vo. In this work, we will employ the RHB+QRPA approach and PC-PK1 effective interaction 
to calculate the half-lives of neutron-rich nuclei, so the value of Vo = 134.0 MeV is replaced 
by 175.0 MeV. Other parameters in Eq. (J3]) take the same values as in Ref. 27]. Therefore, 
the parameters in Eq. (J3]) are not again adjusted with the known half-lives of calculated 
nuclei in this work. 


4 
























In the allowed Gamow-Teller approximation, the /3-decay half-life of an even-even nucleus 
is calculated with 

D 


Tl/2 ~ 


( 4 ) 


where D = 6163.4 ± 3.8 s and qa = 1- The transition strength B m can be directly taken 
from the QRPA calculations. The integrated phase volume /(Z, E m ) is the same as that in 
Ref. ( 27 )]. The /3-decay transition energy E m is calculated with 


Am A n p Aq[{p \ , 


( 5 ) 


where Aq RP a is the QRPA energy with respect to the ground state of the parent nucleus 
and corrected by the difference of the proton and neutron Fermi energies in the parent 
nucleus {t?], 34] and A np is the mass difference between neutron and proton. Since the 
emitted electron energy must be higher than its rest mass m e , the final states must be those 
with excitation energies Aqrpa < A np — m e = A n # (A n # is the mass difference between the 
neutron and the hydrogen atom). 



FIG. 1: (Color online) Nuclear /9-decay half-lives of Fe isotopes, calculated by the RHB+QRPA 
approach with the PC-PK1, compared to the available experimental data. The open squares and 
circles denote the calculations without and with the T = 0 pairing interaction, respectively. 

Fig. ID shows the calculated half-lives of the Fe isotopes in comparison with the available 
experimental data. The results are calculated by RHB+QRPA model with the effective 
interaction PC-PK1. Obviously, the calculated half-lives are very sensitive to the T = 0 
pairing strength, which significantly reduce the calculated /3-decay half-lives. Without the 
inclusion of this pairing channel, the calculated half-lives are systematically longer than the 
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experimental values. With the pairing strength parameter Vo in Eq. (13|) . the half-lives of Fe 
isotopes are reproduced very accurately. 



E (MeV) 

FIG. 2: (Color online) The Gamow-Teller transition probabilities of ,2 Fe calculated without and 
with the T = 0 pairing interaction. The energy threshold A n # for transitions which can contribute 
to nuclear /3 decay is also denoted. 

The nuclear /3-decay half-life is determined by the transition strength as well as the tran¬ 
sition energy which decides the phase volume f(Z , E m ). In order to illustrate the mechanism 
of the influence from T = 0 paring on the /3-decay half-life, the Gamow-Teller transition 
strength distributions of ,2 Fe are shown in Fig. [2] For the calculation without T = 0 pairing, 
the /3-decay half-life of ,2 Fe is mainly determined by the transition at E = —4.95 MeV. This 
transition is dominated by the back spin-flip configuration z/1/5/2 —> 7rl/7/2. Because both 
z/1/5/2 and 7rl/7/2 orbits are partially occupied, the T = 0 pairing could give a remarkable 
contribution to the QRPA matrices related to the z/1/5/2 —)■ 7rl/7/2 pair. When the at¬ 
tractive T = 0 pairing is included, the energy of transition dominated by the configuration 
z/1/5/2 —>■ zrl/7/2 is reduced to E — —7.93 MeV and simultaneously the transition strength 
is also increased. Therefore, the inclusion of T = 0 pairing can reduce the /3-decay half-life. 

For comparison with other theoretical results, nuclear /3-decay half-lives for Fe and Zn 
isotopes are shown in Fig. EJ Our calculations with PC-PK1 well reproduce the experimental 
data of Fe isotopes, including the recent experimental data of ,2 Fe. For Zn isotopes, the 
calculations with PC-PK1 overestimate the half-lives of 78,80 Zn, while it agrees well with the 
experimental half-life of 82 Zn. Since 82 Zn has crossed the shell closure of N = 50, the single¬ 
neutron state vlg 7 /2 is partially occupied. This can lead to a transition dominated by the 
configuration z/lg 7 / 2 —> nlg$/ 2 , which is remarkably influenced by T = 0 pairing. Therefore, 
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FIG. 3: (Color online) Nuclear /3-decay half-lives for Fe and Zn isotopes calculated by RHB+QRPA 
model with the effective interaction PC-PK1. For comparison, the experimental data and theoret 
ical results obtained from the RHB+QRPA calculations with the effective interaction D 3 C* 


the FRDM+QRPA calculations 
DF3a energy density functional 


lj, as well as the continuum QRPA calculations 


m, 


361 ] with the 


371 ] are also shown. 


the half-life of 82 Zn is reduced and then well reproduce the experimental data. For the half- 
lives of ' 8,80 Zn, the influence of T = 0 pairing is much weaker and recent studies have found 
that the tensor force 38] and particle-vibration coupling [39j can improve the overestimation 


of /3-decay half-lives for nuclei with N before the shell closure. For calculations with D 3 C*, 
the results are close to experiment data for 66 ’ 68 ’ 70 p e with a Vo = 125 MeV determined by 
fitting to half-lives of Fe isotopes, while an overestimation of half-lives of 64 ’ 72 Fe is observed. 
For reproducing half-lives of Zn isotopes, Vo is refitted for the calculations with D 3 C*. 
With an enhanced value of Vo = 300 MeV, the half-life of ,8 Zn is well reproduced, while 
the half-lives of 80,82 Zn are still overestimated. Moreover, the results with FRDM+QRPA 
approach systematically overestimate the experimental /3-decay half-lives for both Fe and 
Zn isotopes. This overestimation can be at least partially attributed to the neglect of T = 0 
paring [16|,127]- In addition, by slightly modifying the spin-orbit terms in the DF3 functional, 
the modified functional DF3a has been developed and well describes the experimental data 
on spin-orbit splitting in magic and semimagic nuclei j^]. The available continuum QRPA 
calculations with DF3a agree well with the half-lives of both Fe and Zn isotopes. 

Furthermore, with the isospin-dependent T = 0 pairing strength in Eq. ([3]), we calculate 
the /3-decay half-lives of neutron-rich even-even nuclei with 8 ^ Z ^ 30 by RHB+QRPA 
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FIG. 4: (Color online) Comparison of the calculated nuclear /3-decay half-lives by RHB+QRPA 
model and effective interaction PC-PK1 with experimental data and FRDM+QRPA calculations 
for Z = 8 — 30 even-even nuclei. 


model and effective interaction PC-PK1. Fig. [4] displays the calculated results, experimental 
data, and FRDM+QRPA calculations. It is found that the FRDM+QRPA model system¬ 
atically overestimates nuclear /3-decay half-lives for most of nuclei in this region. However, 
the calculations with RHB+QRPA model generally agree well with the experimental data, 
except for S, Ni, and Zn isotopes, whose half-lives are systematically overestimated. The 
overestimation of half-lives for Ni isotopes can be understood as the same explanation for Zn 
isotopes. When the neutron number of Ni isotope crosses N — 50 shell closure, the half-life 
deviation is remarkably reduced and even reproduce the recent experimental data of 80 Ni [11 ]. 
The overestimation of half-lives for Ni isotopes with neutron number before N = 50 shell 
closure is a common problem in self-consistent relativistic QRPA calculations, which is also 

fin n 

found in RHB+QRPA model with the meson-exchange effective interactions |24|, |25|, [27]. By 
increasing the nucleon effective mass, half-lives of Ni isotopes can be reduced, hence become 
close to the experimental data. However, the increase of m* in relativistic model is limited 
by a realistic description of nuclear matter equation of state and ground-state properties of 
finite nuclei 24), [25]. Therefore, some other effects in mean field framework, such as tensor 
force, and the effects beyond the mean field, such as the particle vibration coupling effect, 
should be investigated for better reproducing nuclear /3-decay half-lives. For S isotopes, it 
was indicated that 40 S and 42 S are deformed 7-soft nuclei, while 44 S exhibits shape mixing 


in the low-lying states 


40 


4l|. Since we ignore deformation in present calculation, the dis¬ 


crepancy between theoretical results and experimental data may originate from the effect of 
deformation for S isotopes. With phenomenological QRPA model, it has been found that 




















the inclusion of deformation indeed reduces the calculated half-lives [ 8 ]. 

In summary, self-consistent proton-neutron quasiparticle random phase approximation 
approach is employed to calculate the nuclear /3-decay half-lives with the relativistic nonlin¬ 
ear point-coupling effective interaction PC-PK1. It is found that the nuclear /3-decay half- 
lives are sensitive to the isoscalar proton-neutron pairing interaction, which can significantly 
reduce the /3-decay half-lives. With an isospin-dependent T = 0 pairing strength, our results 
well reproduce the /3-decay half-lives of neutron-rich even-even nuclei with 8 ^ Z ^ 30 ex¬ 
cept for Ni, Zn, and S isotopes, although the T = 0 pairing strength is not adjusted with the 
half-lives calculated in this work. The underlying reasons for the overestimation of half-lives 
of Ni, Zn, and S isotopes are also discussed in detail. 

This work was partly supported by the National Natural Science Foundation of China 
(Grants No. 11205004, No. 11305161, and No. 11175001), the 211 Project of Anhui Uni¬ 
versity under Grant No. J01001319-J10113190081. 
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